Stretched exponential relaxation and ac universality in disordered dielectrics 
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This paper is concerned with the connection between the properties of dielectric relaxation and ac 
(alternating-current) conduction in disordered dielectrics. The discussion is divided between the clas- 
sical linear-response theory and a self-consistent dynamical modeling. The key issues are, stretched 
exponential character of dielectric relaxation, power-law power spectral density, and anomalous de- 
pendence of ac conduction coefficient on frequency. We propose a self-consistent model of dielectric 
relaxation, in which the relaxations are described by a stretched exponential decay function. Mathe- 
matically, our study refers to the expanding area of fractional calculus and we propose a systematic 
derivation of the fractional relaxation and fractional diffusion equations from the property of ac 
universality. 



I. INTRODUCTION 

Many materials with a disordered structure show a di- 
electric relaxation that is not described by an exponential 
(i.e., Debye-like) decay with a characteristic single decay 
time. Rather the relaxations follow a stretched exponen- 
tial the so-called Kohlrausch- Williams- Watts (KWW) 



</>f,(t) = exp[-{t/ry 



(1) 



with the exponent < < 1 and r a constant. The 
(3 values generally depend on the absolute temperature 
and the chemical composition of the material and typi- 
cally span between 0.3 and 0.8. The relaxation pattern in 
Eq. (1) has been found empirically in various amorphous 
materials as for instance in many polymers and glass- like 
materials near the glass transition temperature (for re- 
view see Refs. H andH, and references therein). Some 
physical models yielding general features of the KWW 
type dielectric relaxation are discussed in Refs. and 
Ifj where one also finds a review of experimental dielec- 
tric relaxation data. Other applications of the KWW 
relaxation function include long-time decay in trapping 
processes^ non-radiative exciton recombination^ and re- 
laxation in sand-piles^ 

The stretched exponential KWW decay function (t) 
can conveniently be considered as a weighted average of 
single-exponential relaxation functions 



(2) 



where the weighting function gp (/j.) is expressible in terms 
of a stable (Levy) distribution^ (Eq. (|26|) below). Be- 
cause of this connection with the statistics of stable laws, 
the KWW relaxation properties can be argued to ap- 
pear naturally through the dynamics, thus challenging 
the common view that the stretched exponential func- 
tion is just a suitable phenomenological fitting tool with- 
out fundamental significance*^ We draw attention to the 



fact that the KWW relaxations arise from a superposi- 
tion of many single-exponential relaxation processes and 
are multi-scale, contrary to a Debye relaxation. 

In this paper we are concerned with the connection be- 
tween the KWW relaxation function in Eq. (1) and the 
anomalous frequency dependence of ac conductivity of 
homogeneously disordered insulators. It has been found 
also partly empirically that insulating and/or poorly con- 
ducting materials with molecular and/or structural disor- 
der exhibit a common conductivity (dielectric) response, 
a phenomenon often referred to as "ac universality" (se e 
review [ToL and original works in Refs. [ll| and Il2t ). 
This response is characterized by low-frequency conduc- 
tivity with very weak or no frequency dependence, and 
a higher-frequency counterpart that follows an approxi- 
mate power law 



a (u) oc a/' 



(3) 



with u'(uj) the real part of the frequency-dependent com- 
plex conductivity <j{uj) and the exponent r\ ranging be- 
tween and 1, and most often between 0.6 and 1, de- 
pending on the material and the absolute temperature. 
The defining feature of ac universality is independence of 
the microscopic details of the disorder and of the nature 
of the charge conduction mechanism operating in the sys- 
tem (classical barrier crossing for ions and/or quantum 
mechanical tunneling for electrons). Note that the power- 
law in Eq. ([3]) is well-defined for frequencies higher than 
the dielectric loss peak frequency. The system-specific 
properties are contained in the coefficient in front of c'(w) 
(not shown in the scaling relation) . The signatures of ac 
universality have been found in materials as diverse as ion 
conducting glasses, amorphous and polycrystalline semi- 
conductors, organic-inorganic composites, ion and elec- 
tron conducting polymers, and doped semiconductors at 
helium temperatures (Ref. 10 for a review). The ob- 
served characteristics including the temperature depen- 
dence of the exponent rj could be reproduced in a model* 3 
in which the conductivity is caused by random motion of 
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FIG. 1: Schematic of cycles and dead-ends of a percolating 
fractal structure. 



charged-particles on a fractal latticed The scale- free be- 
havior of the ac conductivity at higher frequencies could 
be explained as arising from the hierarchic structure of 
the fractal (see Fig. 1) in which multiple cycles and 
dead ends produce effective potential traps influencing 
the motion of the charge carriers. This fractal-based ap- 
proach conforms in spirit with the issue of trap-controlled 
conduction and multiple-trapping transport of charges in 
disordered media (Ref. [HI and references therein). 

The purpose of this paper is to explore the connec- 
tion between ac universality and the stretched exponen- 
tial KWW type dielectric relaxation in disordered di- 
electrics. We shall argue that the two phenomena share 
the statistical-mechanical foundation, and a simple rela- 
tion between the exponents (3 and 7/ will be derived: 



local response of the material the polarization field at 
time t at point r can be written as 



P + ri = 1 



(4) 



The paper is organized as follows. In Sec. II, we show 
that a KWW-type dielectric relaxation implies a power- 
law memory response function for short times and, re- 
lated to this feature, high-frequency dependence of ac 
conduction coefficient with the exponent 77 = 1 — /3. 
We then discuss the route to the power spectrum and 
the occurrence of power-laws for high frequencies in the 
power spectral density (PSD). As a consistency check, 
we re-derive the power spectrum from a Pareto-Levy dis- 
tribution of relaxation times and demonstrate that the 
power-law spectrum is directly related to a power-law 
distribution of relaxation times for the short time scales. 
In Sec. Ill, we present a self-consistent model of di- 
electric relaxation, in which the polarization and elec- 
tric source fields are self-consistently generated by the 
residual polarization-charge density. We show that, if 
the dielectric observes ac universality, the relaxations are 
stretched-exponential for short times. Mathematically, 
our study refers to the expanding area of fractional cal- 
culus and we propose a systematic derivation of the frac- 
tional relaxation and fractional diffusion equations from 
the property of ac universality. We summarize our con- 
clusions in Sec. IV. 



II. RESPONSE FUNCTIONS, SCALING, AND 
POWER SPECTRAL DENSITY 

A. Response functions 

Let a homogeneous, isotropic dielectric be exposed to 
the external polarizing electric field E = E(t, r), which 
we consider as a function of time t and the space co- 
ordinate, r. By homogeneous and isotropic we refer to 
spatial scales larger than the typical scales of the molecu- 
lar/structural disorder. Assuming a linear and spatially 



P(t, r) 



+ OO 



x{t-t')E(t',r)dt' 



(5) 



where x(t~t') is a response or memory function. Causal- 
ity requires that x(i — t') = for t < t'. By consider- 
ing the source field of the form E(t, r) = E(r)S(t) with 
d(t) the Dirac delta function we have the polarization re- 
sponse P(i,r) = E(r)x(i), hence x{t) is the response to 
a delta-pulse in the source field. 

In a basic theory of the dielectric relaxation one is in- 
terested in the polarization response to a field which is 
steady for t < and, then, is suddenly removed at time 
t = 0. One introduces the relaxation function 0(i) as 
the magnitude of the polarization response on an elec- 
tric field with time history E(t, r) = E(r)6*(— t) exp (vt) 
in the limit v — > +0, where 8(t) is the Heaviside step 
function. The infinitely slow exponential growth exp (vt) 
is included in order to satisfy E(£,r) — for t — > —00, 
required by causality. From Eq. §5§ one can see that the 
polarization response is P(t,r) = E(r)0(rj), where 



<j>(t) = 8(-t)(f>(0) exp (vt) + 8(t) U(0) - J X (t')dt' 
and we defined 



0(0) = / x(t'W 

Jo 



(6) 
(7) 



In the literature one commonly writes the value of 0(0) in 
terms of the static permittivity of the dielectric medium 
e(0), i.e., 0(0) = [e(0) - 1] /4vr. The definition in Eqs. © 
and |[7J| is practical for exploring the connection between 
the memory kernel and the relaxation function. For these 
practical reasons, we shall set 0(0) = 1 in the calculation 
below. Performing d/dt on Eq. ^ one obtains 



(8) 



If we postulate the relaxation function in the KWW form 
0(t) = 9(-t) exp (vt) + 6(t) exp [-[t/rf] (9) 
from Eq. §8§ we find 

X (t) = ^- 1 e(t)eM-(t/rf] (10) 

This is a very important expression, because it reveals 
that the stretched exponential relaxation function 0(t) 
leads to a response function x(rj) which behaves like 
a power law t' 3-1 on short time scales t < r, with a 
stretched exponential cut-off on long time-scales t > t. 
This cut-off is necessary for the existence of a finite re- 
laxation function as a response to the step-function elec- 
tric field E(t, r) = E(r)0(— t) exp (vt), since without it 
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the integrals in Eqs. ([5]) and (O will diverge. How- 
ever, if the electric field is an oscillating function or a 
white noise, the response given by Eq. ([5]) will converge 
even with %(£) ~ t^~ x without such a cut-off. The re- 
sponse corresponding to such an "unscreened" power-law 
response function belongs to the class of fractional Brow- 
nian functions^ and will be discussed in some detail in 
Sec. Ill CI By Fourier transforming Eq. ^ we have 

/0 r + OG 

exy{vt)e lu)t dt + / exp[-{t/Tf]e luJt dt 
-oo ^0 

yielding, for v —>■ +0 and uj higher than a non-vanishing 
lower bound (physically corresponding to the dielectric 
loss peak frequency): 



(w) = tQ(wt) + it V(wr) 



i 



(12) 



where Q and V are the Levy definite integrals: 

/■+oo 

Q(z) = / exp (— m' 3 ) cos (uz)du (13) 



r+oo 

V(z) = / exp (— u 13 ) sin (uz)du 
Jo 



(14) 



here expressed as functions of dimensionless frequency 



B. Leading-term approximation 

Series expansion of the Levy functions Q(z) and V(z) 
has been derived and discussed in the literature 
Here we utilize an expansion good for higher frequencies, 
which goes in inverse powers of z: 



Q(z) 



1 r(n/3+l) . nPir 



V(-l)"" 1 

I z np+i r(n+l) 



sin - 



v(z) = J2(-iy 



T(nf3 + 1) nf3n 
■ cos — — 



z np+i r(n + l) 



(15) 



(16) 



From Eqs. (fT5| and (fit))) one can see that the expansion of 
Q(z) starts from a term which is proportional to 
and so does the expansion of V(z) — 1/z. Hence, up to 
higher order terms, 4>(u>) oc uj~( 1+ P\ 

If we now define the frequency-dependent complex sus- 
ceptibility of the dielectric, x(u>) as a Fourier pair with 
X(t), from Eqs. © and (Q2J) we find 

X(w) = iu)<f)(bj) = 1 — lotV(ujt) + io;rQ(wr) (17) 

with the leading term x(u>) oc cj _/3 . With use of the 
Kramers-Kronig relation x(w) oc P J duj'a(uj')/uj'(uj' —uj) 
the scaling of the frequency-dependent complex conduc- 
tivity of the material can be evaluated as 

a{uo) oc uj 1 - (18) 

which reproduces the phenomenological form in Eq. (3) 
with rj = 1 — f3, whence Eq. (5) is obtained. 



FIG. 2: The square amplitude of the frequency-dependent 
complex susceptibility vs the normalized frequency, ujt, for 
different values of the relaxation exponent (3. The plots show 
a transition to a power-law regime for ujt > 1. 



C. Power spectral density 

By Fourier transforming Eq. |(SJ) we find the following 
expression for the PSD of the polarization field: 



S( W )H|P(^r)| 2 )HxM| 2 (|E( W ,r)| 2 > 



(19) 



where the angle brackets (•) denote an ensemble average. 
From Eq. (fl"9|) one can see that |x(^)| 2 is the PSD of 
the polarization field when the driving electric field is an 
uncorrelated white noise signal. In Fig. [2] we have plotted 
(in log- log axes) the |x(^)| 2 function from Eq. (fT7|) . which 
shows a transition to a power-law regime for lot > 1. 

The origin of the power-law regime at the high fre- 
quencies can be found in the short-time behavior of x(i) 
which diverges as i' 3-1 when t — * 0. In fact, the Fourier 
integral 



^exp 



-{t/Tf +iujt]dt (20) 



converges for t/r — > if (3 > 0. This integral is essentially 
the Fourier transform of the power-law function tP~ x with 
a stretched exponential cut-off starting to take effect in 
the integrand for time scales t ~ r or longer (i.e., the 
frequencies uj ~ 27r/r or smaller). For lot sufficiently 
large the integral has already nearly converged before 
the cut-off region in the integrand is reached, enabling 
one to evaluate the susceptibility function as 



ft_ 



t 13 - 1 exp (iujt) dt oc uj^ 13 (21) 



The scaling in Eq. (|21jl complies with the leading term in 
the series expansion of the complex function 1 — zV(z) + 
izQ(z). 

From Eq. (I19p one obtains, for the correspondingly 
high frequencies, S(oj) oc lu~ 213 provided that the input 
driving field is an uncorrelated white noise signal. 

The short-time power- law behavior \(t) ~ t^~ x corre- 
sponds to a memory kernel which defines, by means of a 
convolution with Gaussian white noise, the family of frac- 
tional Brownian functions^ i.e., Xn(t) ~ t H ~ x / 2 where 
the exponent H is often referred to as the Hurst expo- 
nent. The connection with the KWW relaxation function 
exponent is given by 



H = (3-1/2 



(22) 



Depending on the value of H, one distinguishes between 
fractional Brownian motions (0 < H < 1) and fractional 
Brownian noises (H < 0) . The difference is that a noise- 
like function has a stationary quality in the sense its 
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variance does not asymptotically grow with time. From 
Eq. (|2"2"|) one can see that on time scales shorter than r 
the response signal on a white noise source field is a per- 
sistent fractional Brownian noise for < (3 < 1/2 (since 
— 1/2 < H < 0) and an antipersistent fractional Brown- 
ian motion for 1/2 < < 1 (since < H < 1/2). In the 
Debye limit of (3 — > 1, the response signal is an ordinary 
Brownian motion with H — 1/2. Note that (3 = 1/2 cor- 
responding to H = is a special value which separates 
motion-like and noise-like response processes. 

On time scales longer than r, the memory function 
x(t) is strongly influenced by the stretched exponential 
decaying factor. For f3 < 1 the long stretched exponen- 
tial range leads to a deviation from the flat white noise 
spectrum for the low frequencies uj <C 2tt/t and more so 
for smaller (3. This implies correlations on time scales 
( » t (though not yet long-range correlations^ in the 
strict sense, since the latter would require a spectrum 
which diverges faster than logarithmic as lu — > 0). 

From Eq. we have <fi(u>) — —ix{u)/<jj. Hence, the 
power spectrum of the response to a step function driving 
field is 



SH = |0M| 2 = ^- 2 |xH| 2 cx W 



-2(1+0) 



(23) 



The extra factor lo~ 2 in the PSD compared to the PSD 
with an uncorrelated white noise driving field is due to 
the non-stationarity of the step function source signal. 



fi <^ t the parameter r / fi corresponds to the tail of the 
Levy distribution, which is approximated by the Pareto 
inverse-power distribution,— i.e., 



-•13,-1 



(r/ / i)«^(r/ / i)-( 1 +« 



(27) 



where Ap is a normalization parameter and < (3 < 1. 
The distribution of relaxation times is then a pure power 
law: 



gpifJ-) oc (i 2 /i 1+/3 oc fi ^ 



(28) 



Now let us demonstrate directly the connection between 
the power-law distributed relaxation times for \i <C t 
and the power-law power spectra for u> 3> 2tt/t. By 
considering to 3> 2t:/t we can define a cross-over time 
scale a ~ Q/u> with a constant much larger than 1. 
We can then split the integration in Eq. ((25)) into a tail 
integral J Q dfi and a core integral J a dfi. In the core 
integral we have ufi > O ^> 1 and we can neglect 1 in 
1 — iufi, then integrate the distribution function gpifi) 
through all fi to find, with J gp(fi)dfi — 1, 



a 1 - UJ 



(29) 



Thus the core integral compensates for — i/ui in Eq. (|25[) . 
leading to 



D. The physical origin of power- law spectra and 
stretched exponential relaxation 

It is instructive to re-derive the power spectrum in 
Eq. (f2"3")l from the distribution of relaxation times as char- 
acterized by Eq. (2), then from the weighted superposi- 
tion of the Debye single-exponential relaxation processes. 
From Eq. (|11[) we have 

■ poo poo 

<f>(u) = +/ e iut dt / Q (n)exp(-t/n)dn (24) 

w Jo Jo 

where we expanded the stretched exponential over the 
partial relaxation times, /i. Changing order of integra- 
tion in Eq. (|24p and Fourier transforming the exponential 
function exp(—t/fi) one finds 



H = — 

U) 



fl 



1 — iuifi 



(25) 



The weighting function Qp(fJ>) is given by Eqs. (51d) and 
(55) of Ref. [5j where one replaces the exponent a with 
(3, the time constant T with r, and the variable fj, with 
r/fi. In our notations: 



0/3 (/«) 



-•13, 



(26) 



where L^_i is the Levy distribution function with skew- 
ness parameter —1 (Ref. 19). For short relaxation times 



1 



-gp(/i)dfj, 



(30) 



In the tail integral we have r/fj, S> 1, and we can sub- 
stitute Qff(p) from Eq. (|26p and utilize the power-law 
distribution in Eq. (|28|) to obtain 



( W ) oc f -JL—dpL w 
Jo 1 - IWJi 



(31) 



with ^ = w/x varying from to Q. The integral on the 
right of Eq. (f3"Tj) converges to a constant, while the scaling 
factor in front yields the PSD in Eq. ([25]). 

Our findings so far can be summarized as follows: The 
stretched exponential relaxation corresponds to a power- 
law memory response function for short time scales t < r. 
For a white noise driving electric field this in turn im- 
plies a fractional Brownian polarization response signal 
on these time scales. These features are intimately con- 
nected with the description of the stretched exponential 
relaxation as a superposition of exponentially decaying 
signals with a distribution of characteristic relaxation 
times fi which are power-law distributed for /i <C t. 

Power-law distribution of durations of relaxation 
events as responses to external perturbation is a hall- 
mark of systems in states of self-organized criticality. 20 
The discussion presented in this section supports the hy- 
pothesis that dielectrics exhibiting stretched exponential 
relaxations are in a self-organized critical state. 
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E. Fractional-derivative representation 

If one wishes to calculate the spectrum S(u>) from the 
dynamics of the relaxation process, the procedure is to 
replace the memory function in Eq. |(SJ| by its KWW 
representative in Eq. (fT0|) to give 

^^jyj^^-^r (32 , 

The integration in Eq. (f32|) is singular at t' — t, but the 
integral converges if < j3 < 1. As mentioned in Sec. Ill Al 
the stretched exponential cut-off factor in the integrand 
is necessary for convergence of the integral if E(t', r) is 
a step function. However, for a white noise electric field 
the integral converges even if this cut-off is removed by 
replacing the exponential factor with unity. The effect 
this removal will have on the nature of the response sig- 
nal is to introduce long-range statistical dependence on 
time scales longer than r, while the statistical properties 
on scales shorter than r are unaffected. Focusing on the 
shorter time scales we shall proceed with the cut-off fac- 
tor removed. By applying d/dt to both sides of Eq. (j3"!?|) 
we have 

T^P(t,r) = _ oc D ( 1 - /3 E(t,r) (33) 

with Tp = r' 3 /r(l + 0) a constant and _ 00 Dj - ' 3 the 
so-called Riesz fractional derivative^! which is defined 
through 

-ai-'/M-fWs/L*^^ (34 » 

where f(t,r) is differintegrable for t' — > t at point r. 
Along with the Riemann-Liouville derivative (to be dis- 
cussed in some detail below) the Riesz derivative in 
Eq. offers a fractional generalization of ordinary 

derivative d/dt, which can be thought as a special case of 
-ooDj - ^ with the integer (3 — and 1-/3=1. A Fourier 
transformed Riesz derivative _ 00 Dj^' 5 is (— icj) 1- ' 3 , in 
analogy to the transform of d/dt. Performing a Fourier 
transform of Eq. (j3"3"|) we get 

- ttjP(cj, r) oc (-iw) 1 "' 3 E(w, r) (35) 

where we omitted T p for simplicity. From Eq. (|35p one 
immediately recovers the PSD of the polarization field. 
In the integrated form Eq. (|33|) reads 

T P(t,r) =_ 00 D t - /3 E(t,r) (36) 

with „ 00 D 4 "' 3 acting as 

The operator in Eq. (|3~T|) is known as the Riesz fractional 
integral. With (3 = H +1/2 and E(t, r) a Gaussian white 



noise the integration in Eq. p6|) generates a fractional 
Brownian function (i.e., antipersistent fractional Brow- 
nian motion for 1/2 < j3 < 1 and persistent fractional 
Brownian noise for < (3 < 1/2). The origin of frac- 
tional Brownian type polarization response to a white 
noise electric source field can be imagined as arising from 
non-random motion of dielectric molecules driven by the 
uncorrelated external forcing. 

III. FRACTIONAL KINETIC EQUATIONS 
A. Self-consistent dynamic-relaxation model 

We have assumed in Sec. II that the electric field 
E(i',r) has the external origin and we have investigated 
the properties of the polarization response field for dif- 
ferent forms of the electric source field (i.e., the white 
noise and step function driving fields). In this section 
the assumption of the external origin will be relaxed and 
we shall consider E(i', r) as the inherent field of the po- 
larization charges. 

We propose that if, in a dielectric medium, the prop- 
erty of ac universality is verified, then the decay of polar- 
ization goes as a KWW stretched exponential relaxation 
function. We formulate the problem as an initial-value 
problem for time non-local polarization field of the form 

r+oo 

P(t,r) = P(0,r)+ / X (t-t')V{t',v)dt' (38) 
Jo 

where P(0, r) is the initial polarization. We shall assume 
that the response processes to all the external driving 
fields have been accomplished by time t — and we are 
interested in the self-consistent dynamics of relaxation 
when the residual polarization field is essentially the re- 
sponse to the residual electric field due to the polarization 
charges. 

Let p(t, r) be the density of the polarization charges 
at time t at point r. The function p(t, r) is defined as 
the mean density of the electric charges in a physically 
small volume around r such that the highly fluctuating 
molecular densities are averaged out. One assumes that 
there is a length scale separation between the microscopic 
molecular scales and the length scales on which the mean 
density varies. On length scales comparable to or shorter 
than the molecular scales, the dynamics of relaxation 
must be described statistically. We shall return to this 
issue in Sec. Ill C. 

Assume for simplicity that there are no external 
charges, i.e., the total charge of the medium is zero. The 
polarization charges are then the only source for the po- 
larization and electric fields, i.e., 

V • E(i, r) = Anp(t, r) (39) 

and 

V-P(i,r)=-p(f,r) (40) 
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Hence 



V-D(t,r) 







(41) 



where D = E + 4irP is the electric displacement in the 
medium. The density of the polarization/relaxation cur- 
rents is defined as the time derivative of the polarization 
field, i.e., j = dP/dt. From Eq. (38]) we have 



j(t,r) = / a(t-t')B{t',r)dt' (42) 



where 

a(t-t') = ^ X (t-t') (43) 

is a new memory function. A Fourier transformed a(t) 
defines the frequency-dependent complex conductivity 
of the medium, which is related to the susceptibility 
function via u(lu) = —iu>x(^>)- Causality requires that 
a(t — t') = for t < t' so that the integration in Eq. |g2]) 
is non-trivial only in the window < t' < t. Performing 
d/dt on Eq. ([40]) we obtain the continuity equation 



0_ 



p(i,r)+V.j(t,r) =0 



(44) 



with j(£, r) in Eq. (|42|) . In the following we proceed with 
the derivation of the KWW stretched exponential func- 
tion. 



B. KWW decay function and fractional relaxation 
equation 



By Laplace transforming Eqs. ([39]). (l42|) . and (|44j) we 



get 



V • E(s,r) = 4?rp(s,r) 



j(s,r) = cr(s)E(s,r) 



(45) 



(46) 



sp(s,r)-p(0,r)+V-j(s,r)=G (47) 

If we apply V- on Eq. (|46|) and utilize Eq. (f45|) we find 

V-j(s,r)=<T(s)V-E(s,r)=47ra(s)p(s,r) (48) 

When substituted into Eq. (j47|) this yields 

sp(s,r) + 4vrcr(s) / 9(s, r) = p(0, r) (49) 

We assume that the sample observes ac universality, i.e., 
<t(s) = as v with a a constant and rj a fraction between 
and 1. This power-law form is just the Laplace version 
of Eq. ([3]). We shall use this form as an approximation 
of <r(s) for the corresponding (short) time scales, t ~ 1/s. 
Equation (l4"9"[) becomes 



sp(s, r) + 47ras' 7 p(s, r) = p(0, r) 



(50) 



Separating variables we write p(s, r) = (f>(s)ip(r) with (p(s) 
the Laplace transform of the relaxation function 4>(t). 
From Eq. ([50]) it follows that 



(s) = 



1 



s + Airasv 



(51) 



where the initial condition </>(0) = 1 has been applied. In 
the time domain, 



^ ^ 2-ki J_ loo 



S + T 



-/3 s l-/3 



da 



(52) 



where we used the notations (3=1 — 7] and r ° = 47ra. 
Equation l|52p is exactly the definition of the Mittag- 
Leffler function [-(t/r) 73 ] (Eq. (B.l) in Appendix 
B of Ref. l2~ih . The Mittag-Leffler function has the series 
expansion 



E r 



-(t/rY 



n=0 v f 



,n8 



1 



(53) 



For short times, this expansion goes as a stretched expo- 
nential, i.e., 



E/j «exp 



r(/3 + i) 



(54) 



This closed analytic form replicates the KWW relaxation 
function in Eq. (1). 

We now derive a dynamical relaxation equation from 
the dispersion relation (|50f . The key step is to notice that 
the power-law s n with < 77 < 1 is the Laplace trans- 
formed Riemann-Liouville derivative,— which is defined 
through 



oD?/(t,r) 



1 



d 



T(l - rj) dt 



(t - t'Y 



(55) 



The Riemann-Liouville derivative differs from the Riesz 
derivative in that the integration starts from t = and 
not from t — — 00. Despite some particularities^ of the 
composition rules and initial-value terms both vehicles 
share the property of being well-defined fractional ex- 
tensions of the ordinary differentiation. In the limit of 
7] — > 1, the fractional derivative in Eq. (J5SJ) reduces to 
ordinary time derivative, d/dt. Replacing s n by oD^ in 
Eq. ([50)1 we write, with t~@ = Aira, 



d_ 

dt 



p{t,T) 



-r-" D t 1 - / y*,r 



(56) 



It is instructive to derive Eq. (|56|) directly from Eq 
Applying V- to Eq. ([38]) and utilizing Eqs. (J39J) and (gDJ) 
we find, with p(0, r) = - V • P(0, r), 

p(t, r) = p(0, r) - 4tt / X (t - t')p(t', r)dt' (57) 



which is a closed integral equation for p(t, r). Note that 
the causality condition — i') = for t < t' has allowed 
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us to set the upper limit of integration to t. By time 
differentiating Eq. ([57]) we get 



0_ 

Of 



p(t, r) = -4tt^ J X (t - t')p(t', r)dt' (58) 



The memory kernel is calculated from Eqs. © 

and (|52|) . Here we are interested in the short-time be- 
havior, which can be most readily evaluated from the 
power expansion in Eq. (|53|) to give 



x(t) 



(59) 



Of course, the scaling in Eq. (|59[) may be obtained as the 
inverse Laplace transform of the complex susceptibility, 
X(s) cx c(s)/s oc s - ^. When the power-law in Eq. (|59|) is 
substituted to Eq. (|58[). the Riemann-Liuville derivative 
Dj~ is built, whence the dynamical Eq. (jB"6")) follows. 

Equation (IB6|) with the Riemann-Liouville fractional 
derivative is the canonical form of the fractional relax- 
ation equation^ The short-time behavior of the solu- 
tion of the fractional relaxation equation is the KWW 
stretched-exponential decay function, in accordance with 
the power expansion of the Mittag-LefRer function. 

Our main finding in Sec. Ill B is: The KWW stretched 
exponential decay function can be derived from the ba- 
sic electrostatic equations under the additional assump- 
tion that ac conduction coefficient behaves as a fractional 
power of frequency. 



is defined as the frequency-dependent complex diffusion 
coefficient, T>(w). The value of T>(u) can be expressed in 
terms of the ac conduction coefficient as 



T 



cr(w) 



(61) 



where e denotes the carrier charge, n their number den- 
sity, and T the absolute temperature. In the zero- 
frequency limit, Eq. (|6ip reduces to the conventional 
Einstein relation between the diffusion constant and the 
dc conductivity. Based on Eq. ^ we can argue that, 
2?(u>) oc u) 71 for the corresponding high frequencies. 

We now turn to demonstrate that, if the sample 
observes ac universality, the microscopic dynamics of 
charges is described by a fractional extension of the dif- 
fusion equation. 

By Laplace transforming Eq. (fST))) we get 



j(s,r) = -D(s)Vp(s,r) 



(62) 



where T>(s) is the Laplace transform of T>(t). When sub- 
stituted into the continuity Eq. (|4"T)) this yields 



(63) 



sp(s,r)-p(0,r) =©(s)V 2 p(s,r) 



Adhering to the power-law form <r(s) = as 17 from Eq. (|6"Tj) 
we have T>(s) = As n with A = aT/ne 2 . Utilizing this in 
Eq. ([63]) we write 



sp(s,r)-p(0,r) = s"V 2 p(s,r) 



(64) 



Fractional diffusion equation describing 
sub-diffusion 



Our purpose now is to contrast the electrostatic de- 
scription of the decay of polarization with a statistical- 
mechanical description of the dynamics of charged parti- 
cles on the microscopic scales of the molecular motions. 

Microscopically, the relaxations are due to the motion 
of charges which interact with the fluctuating molecular 
environment. As a model approximation, we shall rely 
on the hypothesis of trap-controlled conduction and dif- 
fusion, in which the transport occurs as a result of hop- 
ping of charged-particles between the localized states. If 
the hopping has a characteristic time, then the transport 
is described by a Markovian chain process. In a more 
general situation there is a distribution of waiting or res- 
idence times between the consecutive steps of the motion 
and the Markovian property is invalidated. The system 
response to a charge-density perturbation is then a flow 
with memory: 



3(t,r) 



V(t-t')Vp{t',r)dt' 



(60) 



which goes against the concentration gradient as due to 
Fick's law. Here, V(t — t') is a memory function, such 
that V(t -t') = for t < t' . A Fourier transformed V(t) 



where we set A = 1 for simplicity. In the time domain, 
Eq. (f64|) reads, with 77 = 1 — /?, 



^p(t ) r) = D t 1 -^V 2 p(* ) r) 



(65) 



Equation (f6"S"|) is the canonical form of the fractional 
diffusion equation describing sub-diffusion^ with (3 
the fractal dimension in time. 22 In various settings, 
this equation has been derived and discussed in the 
literatureM* 2 ^* 2 ^ 

The characteristic function of the fractional diffusion 
Eq. (|65p obeys the fractional relaxation equation 



aP(*,k) 



-k^D-i-^k) 



(66) 



where k is the wave- vector in the ambient real-space. 
From Eq. (|54|) one can see that the short-time behavior 
of p{t, k) is the stretched exponential 



p(t, k) w exp 



k 2 ^ 



r(/3 + i) 



(67) 



With this observation we conclude the analysis of connec- 
tion between the KWW stretched exponential relaxation 
function and ac universality in disordered solids. 
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IV. SUMMARY 

We have analyzed the properties of dielectric relax- 
ation and ac (alternating-current) conduction in disor- 
dered dielectrics. We have discussed the route to the 
statistical mechanical foundation and suggested physical 
models that might connect to the typically observed dy- 
namical characteristics. We proposed a self-consistent 
model of dielectric relaxation, in which the polarization 
and electric source fields are self-consistently generated 
by the residual polarization-charge density. This self- 
consistent approach has led us to a systematic derivation 



of the fractional relaxation and fractional diffusion equa- 
tions from the property of ac universality. Our results 
support the hypothesis that dielectrics exhibiting ac uni- 
versality and stretched exponential relaxations are in a 
self-organized critical state. 
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